Long-range potential fluctuations and 1/ / noise in hydrogenated amorphous silicon 



m 
O 
O 

(N 
-i— > 
O 

O 
o 

(N 



S3 

S3 
i 

C/3 



I 

S3 
O 

o 



(N 
> 
O 
oo 

o 

(N 

o 

c3 



S3 
O 

o 



X 
S3 



B. V. Fine 1 - 3 *, J. P. R. Bakker 2 , and J. I. Dijkhuis 2 

1 Spinoza Institute and 2 Debye Institute, Utrecht University, P.O. Box 80000, 3508 TA Utrecht, Netherlands 
3 Max Plank Institute for the Physics of Complex Systems, Noethnitzer Str. 38, D-01187, Dresden, Germany 

(Dated: June 4, 2003) 

We present a microscopic theory of the low-frequency voltage noise (known as "1//" noise) in 
^tm-thick films of hydrogenated amorphous silicon. This theory traces the noise back to the long- 
range fluctuations of the Coulomb potential created by deep defects, thereby predicting the absolute 
noise intensity as a function of the distribution of defect activation energies. The predictions of this 
theory are in very good agreement with our own experiments in terms of both the absolute intensity 
and the temperature dependence of the noise spectra. 



PACS numbers: 71.55.Jv, 72.70. 



I. INTRODUCTION 



73.50.Td, 73.61. Jc 



When electric current flows through a resistor, strong 
low-frequency voltage noise superposed on the thermal 
Johnson-Nyquist noised is generally observed. The 
spectrum of that excess noise has shape close to 1//, 
where / is frequency 3 -. In recent decades, the model, ac- 
cording to which the 1// noise is produced by an ensem- 
ble of two-state systems having broadly distributed acti- 
vation energies (BDAE)!*^, has enjoyed a lot of success. 
However, the BDAE model addresses neither the origin 
of the two-state systems nor the noise mechanism, thus 
leaving the absolute noise intensity as an adjustable pa- 
rameter. The lack of direct microscopic calculations pre- 
dicting the absolute noise intensity, while being a glaring 
theoretical gap as such, also fuels an old but still contin- 
uing debate over the question of what actually fluctuates 
— the number of carriers or their mobility 3 . 

In this work, we present a theoretical and experimen- 
tal investigation of 1// noise in hydrogenated amorphous 
silicon (a-Si:H). Our theory links the noise to the fluc- 
tuations of the number of carriers, predicts the absolute 
noise intensity, and also allows us to extract novel and 
useful information about the defects in this important 
material^. 

This theory can be outlined as follows: Hydrogenated 
amorphous silicon has a significant number of deep de- 
fects known as dangling bonds. Thermal fluctuations of 
the numbers of electrons occupying the dangling bonds 
cause long-range potential fluctuations, which give rise 
to fluctuations of the local densities of carriers, which, in 
turn, lead to resistance fluctuations, which, in the pres- 
ence of current, manifest themselves as voltage noise. It 
is the inclusion of the long-range potential fluctuations 
into the above scheme that distinguishes our treatment 
from many similar theoretical proposals. 



II. EXPERIMENTAL SETUP 

The 1// noise in a-Si:H has been studied in the lit- 
erature in a variety of experimental settings (see e.g. 
RefsJi 8 ! 9 ! 10 ! 11 ! 12 ! 13 ! 14 ! 15 ), and exhibited certain features, 



which depend on numerous details of each particular ex- 
periment. In this work we present a fully developed quan- 
titative study of only one situation, which corresponds to 
our own experiments. Other experimental settings will 
only be discussed briefly in the end of the paper. 

The present study is focused on an n — i — n film of 
a-Si:H, where n stands for a 40nm-thick electron doped 
layer, and i for an undoped layer of thickness d = 0.91/im. 
The n — i — n structure is grown by plasma enhanced 
chemical vapor deposition (PECVD) on a highly con- 
ductive wafer of crystaline silicon. The contact layer on 
the top of the structure consists of a 30nm-thick film of 
titanium followed by 30nm-thick film of copper. The film 
has area A = 0.56cm 2 . It is thermally annealed and then 
protected from light. We observe and analyse the voltage 
noise spectra at frequencies / = 1 -f- 10 4 Hz and tempera- 
tures T — 340 -7- 434 K in the presence of electric current 
flowing perpendicular to the plane of the film. Other 
details of our experimental setup are described in Ref. 1 - . 



III. FORMULATION OF THE THEORETICAL 
PROBLEM 

Now we turn to the theoretical derivation of the noise 
spectrum under the above experimental conditions. The 
central and, presumably, quite general part of this deriva- 
tion is Section llVl Most of the rest is specific to a-Si:H 
and to the experimental setting considered. 

Our goal is to compute voltage noise spectrum Sv{f) 
expressed as: 



Sv(f) 
V 2 



= 4 



C v (t)cos{2irft) dt, 



where V is the applied voltage, and 

(SR(t)5R(0)) 



C v {t) 



R 2 



(1) 



(2) 



Here, Rq is the average resistance of the film, and 6R(t) 
is the equilibrium resistance noise. As usual 3 -, the link 
between the resistance noise and the voltage noise is es- 
tablished experimentally by observing that Sy(f) oc V 2 . 
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FIG. 1: (a) Sketch of the density of states in the center of 
the undoped layer, (b) Band bending profile. The gray stripe 
represents the "uniform resistivity layer" defined in the text, 
(c) Cartoon of a deep defect surrounded by medium-range 
structural disorder. Note: e<j>(t, r) fluctuates on a much longer 
lengthscale and with much smaller amplitude. 



The density of states of undoped a-Si:H is shown in 
Fig. QJi. It is characterized by a band gap of 1.8 eV 
between the mobility edges E v and E c in the valence 
and conduction bands, respectively. The proximity of 
the n-layers induces band bending in the undoped layer 
as shown in Fig Around the center of the film 

E c (z) = E c q — f3z 2 , where the z-axis is directed perpen- 
dicular to the film, E c o = p + 0.63 eV, p is the chemical 
potential, and (3 = 1.6eV//xm 2 . The difference E c o — p 
was obtained experimentally as the conductivity activa- 
tion energy, while the value of (3 was found numerically 
(cf.ll). 

Since p is significantly closer to E c than to E v , 
the resistivity p is inversely proportional to the den- 
sity of conduction electrons n e . Keeping in mind that 

n e cx exp ^— E k c ~T ^ , we obtain p(z) = A exp ( r"^ ) ; 
which is strongly peaked around z = 0. Here A is a pro- 
portionality coefficient, and ks is the Boltzmann con- 
stant. The total film re sistan ce can then be found as 

R = \Jp(z)dz^em^r. 

To simplify the treatment we replace the actual profile 
of p(z) by a layer of constant resistivity po = p(0) spread- 
ing along the z-axis between —zt and zt (see Fig. ^p). 
We require the resistance of this constant resistivity layer 
to be equal to Rq, which gives 



the screened Coulomb potential (f>(t, r) created by deep 
defects: 



<^,r) = > exp 



r - a. 



(4) 



Here Aq^i) is the fluctuation of the ith defect charge 
with respect to its average value, the position the de- 
fect, r s% the screening radius of that defect, and e = 12 
the dielectric constant. The defects causing the potential 
fluctuations may be located outside of the constant resis- 
tivity layer. [Note: the above formula does not include 
large static component of the random Coulomb poten- 
tial.] 

When the local potential cf>(t, r) fluctuates, the mobil- 
ity edge tracks it, i.e. E c (t, r) = E c o + ecj)(t, r), where e is 
the electron charge. Since p does not shift with e<fi(t,r), 
the density of conduction electrons re-equilibrates follow- 
ing E c (t,r) on the timescale of electron drift from the 
uniform resistivity layer to the n-layers (and then to the 
contact layers). Because of the strong band bending, the 
drift takes less than 10 _7 s, i.e. the re-equilibration is 
effectively instantaneous on the timescales of the noise 
studied (y- ~ 1CP 4 -j- 1 s). The resistivity can thus be 



rewritten as p(t, r) = X exp ^ 



£c(t,r)- 
k B T 



Assuming for 

a moment [and proving later] that | e</!>(t, r) | <C ksT , we 
expand p(t, r) = p + Sp(t, r), where p = A exp , 
and 



k R l 



(5) 



For Sp <C po, the fluctuation of the total resistance is 
SR(t) = Jjp(t,r)d 3 r, 



(6) 



where V is the space inside the constant resistivity layer. 
Substituting R = and SR(t) given by Eq.® into 

Eq.© and then using Eq.JSJ), we obtain 



C v (t) 



2k B Tz T A 



d 3 r / d 3 r'((/>(i,r)0(O,r')). 
v Jv 



(7) 



DESCRIPTION OF DEFECTS 



(typical value: 0.13/im). 




(3) 



IV. CORRELATION FUNCTION 

The resistivity within the "constant resistivity layer" 
still fluctuates as a consequence of the fluctuations of 



In order to evaluate ((/>(t, r)<fi(0, r')), we have to de- 
scribe the deep defects in the undoped layer of a-Si:H 18 . 

There exist both theoretical arguments and experi- 
mental evidence^ indicating that the concentration of 
defects in the undoped layer should be greater in the 
vicinity of the n-layers and then decay towards the cen- 
ter on the lengthscale of about 0.5/im. Here, however, 
in order to simplify the theoretical treatment, we as- 
sume that the defect concentration across the entire un- 
doped layer has a constant value, which we estimate as 18 
hd = 6x 10 15 cm" 3 . 



3 



Each defect has 4 possible states: one zero-electron 
state D + , two one-electron states D°, and one two- 
electron state, D~. We assume a Gaussian probability 
distribution 



p(E) 



1 



i^AE 



(8) 



for energy of an electron occupying the D° state. Here 
E a = fj, - 0.22 eV and AE = 0.15 eV. Energy E, as such, 
is associated only with the D + <-> D° transition. The 
D° <-* D~ transition, which requires capturing the sec- 
ond electron by the same defect, is characterized by en- 
ergy E' = E + U , where U = 0.2 eV is the correlation en- 
ergy. Therefore, the corresponding probability distribu- 
tion is 



p'(E') 



1 



^AE 



(E'-E -U) 2 



(9) 



We shall treat the D + <-* D° transitions independently 
from the D° <-> D~ transitions, which is justified as long 



as exp 



u 

k B T 



< 1. 



The noise spectrum obtained later will depend only 
weakly on defect parameters no, Eq, U and AE, mainly 
through a weak dependence on one combination of them 
— the density of defect states at the chemical potential. 
That combination, in turn, is not very sensitive to the 
choice of E and U. 

Using the equilibrium description^ of the D + , D° and 
D~ states, we obtain the mean squared charge fluctua- 
tions for the D + <-> D° transition: 



(A^i^+^^sech 2 



E 



2k B T 



where c = | ln2. For the D° <-> D transition: 



(A 9 2 (S'))--7e 2 sech 2 '^' M 



2k B T 



(10) 



(11) 



In order to escape from a deep defect, an electron 
should reach the mobility edge E c . However, the acti- 
vation barriers Eb (indicated in Fig^;) can vary as a 
result of the medium-range disorder of the amorphous 
structure (on a lengthscale of 1 -f- 10 nm). We assume a 
Gaussian probability distribution for the values of Eb- 



P(E B ) 



1 



2-eAEb 



exp 



(Eb — Esq) 2 
2AE% 



(12) 



where Ebo and AEb are to be extracted from the ex- 
perimental spectra. These are the only two adjustable 
parameters in our treatment. They will affect the spec- 
tral shape, but not the integrated noise intensity. 

The fluctuation rate for the D + <-> D° transition is2i 



1 



t+(E,E b ) 



LU (E B ) 



1 _££ 

2 6 " 



(13) 



and for the D° «-> D transition: 
1 



t-(E',E b ) 



= ljo(Eb) 



1 

— P k B T 

2 



(14) 



where the attempt frequency is estimated from the de- 
tailed balance condition as 



uo(Eb) = v t h cr N C {E B ) 



(15) 



Here, v t h = y^StcsT/m* is the thermal veloc- 
ity of conduction electrons; to* their effective mass 
(to* w 0.4to 6 = 3.64 • 10 _28 g); a = 10" 15 cm 2 the cross- 
section of electron capture by D + or D° defect; and 
Af c (E B ) ~ N C (E B ) k B T the "concentration" of ther- 
mally accessible electronic states with energies above 
E B ; N c (E b ) = N c (E c0 )s/(E B - E c0 + e c )/e c the empir- 
ical fit to the density of states above the mobility edge; 
N c (E c0 ) = 4 x 10 21 eV _1 cm- 3 ; and e c = 0.02 eV. The 
typical value of Uo(E B ) is then 10 13 s _1 . 

Next we evaluate the screening radius. Since the 
screening by the conduction electrons in the undoped 
layer can be neglected in view of their very small con- 
centration (10 10 -r 10 13 cm~ 3 ), two other screening mech- 
anisms should be considered, namely: (i) by the n-layers 
together with the contact layers; and (ii) by the defects 
in the undoped layer. 

We describe the first mechanism as a perfect screening 
by metallic surfaces. In other words, we assume that a de- 
fect in the undoped layer is screened by the infinite set of 
its mirror images constructed with respect to the mirror 
planes located at z = ±d/2. Although the screening law 
due to this mechanism is not exponential, we obtain the 
best possible value for r s entering Eq. J3J| as the distance 
from the defect to the point, where the screened potential 
is factor of e = 2.718... smaller than the bare Coulomb 
potential. We have found numerically that r s has angu- 
lar dependence, which is such that, for a defect located at 
2 = 0, r sx = r sy = 0.54d, and r sz = 0.38d. We then ap- 
proximate the screening radius due to this mechanism as 



r s i = (r 



1/3 



0.48d. 



(16) 



The second screening mechanism is similar to that of 
Thomas-Fermi. It can be described as follows: When the 
potential fluctuates due to charge fluctuation on a given 
defect, the energies E and E' of the neighboring defects 
become shifted. In response, those neighboring defects 
change their occupation numbers thus screening the po- 
tential of the defect, which originally caused the fluctu- 
ation. The equilibrium description of this mechanism 22 
results in the screening radius 



ec l \j Aire 2 n D v(Ty 



(17) 



where 



v{t) - ^f{L mdE+ jy^} ; (i8) 



4 



p(E)=p(E) S ech 2 ^^-£-c), (19) 

and 

p'(E') = p'(E') sech 2 + c) . (20) 

The value of v(T) depends on T very weakly. The prod- 
uct n,£)v(0) should be recognized as the density of the 
defect states at the chemical potential. 

The equilibrium description of screening is applicable, 
when the sources of the Coulomb potential are static, and 
thus all the defects can contribute to screening. However, 
if a given defect fluctuates on time scale r, its poten- 
tial can only be screened by other defects fluctuating on 
timescales not slower than t. To take this observation 
into account, we simply multiply n B in formula l|17l) by 
6(r), the fraction of defects that fluctuate fast enough to 
take part in the screening of a defect characterized by the 
fluctuation time r. We estimate that fraction as 

Kr) = 4ldv(T) {le + (r)P( E ) P ( E b) dEdE B 

+ S £ _ {T) p'{E')P(E B )dE'dE B ), (21) 

where £+ (r) is the integration region limited by condition 
t + (E, Eb) < t, and £-(t) the integration region limited 
by t-(E',Eb) < t. The screening radius due to the 
second mechanism thus becomes 

rs2(T) = \/ AneHnl(T)b( T y (22) 
I 



For the slower fraction of defects, r S 2(r) approaches 
r eq ~ 0.18 /xm, which is shorter than r s i = 0.43 /im. At 
the same time, for the faster ones, r S 2 is very large, i.e. 
they are predominantly screened by the first mechanism. 
We, therefore, approximate the combined effect of the 
two mechanisms by using the following expression for the 
screening radius: 

r- 1 (r) = r 7l 1 +r:i(r). (23) 
The typical value of r s given by Ea. (|23J) is r* — 0.2 /im. 



VI. EVALUATION OF THE NOISE SPECTRUM 

Now we are in position to evaluate (4>(t, r)(f>(0, r')) with 
4>(t,r) given by Eq.(gJ. 

Since an electron emitted by one deep defect is most 
likely absorbed not by another defect but by the contact 
layer, it is appropriate to assume that different defects 
fluctuate independently, which implies that, for i =/= j, 
(Aqi(t)Aqj(0)) = 0. Keeping also in mind that 

(A qi (t)A qi (0)) = (Ag 2 ) exp(-t/ri), (24) 

we take the disorder average of <j>(t, r)(p(0, r'), i.e. we 
replace the discrete summation by the integrations over 
the relevant probability distributions. This gives 



(^,r)0(O,r')) = ^ J dE B P{E B ) {j dEp(E) F(r, r', r + (E, E B )) {Aq{Ef)+ exp (- 7T71 b^) 

+ JdE'p'(E')F(r,r',r-(E',E B )) (Ag(£') 2 >- exp (- t _ {e \ Eb) )} ■ (25) 



where 

/ |r— al + lr 7 — a| \ 

J |r- a| |r' - a| 
«2 7 rr s (r)exp(-^). (26) 

The approximation in Ea. (|26|l can be justified by observ- 
ing that the main spatial dependence of the integral for 

F(r, r', r) has form exp ^— ^*~(t) ) ' an< ^ then the natural 
choice for the prefactor is F(0,0,t) = 27rr s (r) (obtained 
by direct integration). 

From Ea. l|25|l . the value of | e<^>(£, r) | can be estimated 
as: ev /(0 2 (O,O)) w e 2 e- 1 ^2Ttn D v{T)k B Tr* s - 3.5meV. 



Since k B T ~ 30meV, the assumption |e0(i,r)| k B T 
made earlier was adequate. 

Finally, we evaluate Cy (t) by substituting Eq. I|25|) into 
Eq. , and then take the Fourier transform Q to obtain 

Sy(f) 27T 2 e 4 n D 
V 2 e 2 (k B T) 2 Az 2 

x J dE B P{E B )^J D{},T+{E,E B ))p{E)dE 

+ / D(f,T^(E',E B ))p'(E')dE'\, (27) 
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where 



D(f,r) 



}(r) 



4z 7 



3r s (r)+e fTPT (2z T + 3r s (r)) 



1 + 4tt 2 / 2 t 2 



,. (28) 

Recalling the approximations associated with (i) the 
assumption of the constant resistivity layer; (ii) the as- 
sumption of the constant profile of defect concentration; 
(iii) the treatment of the screening radius; and (iv) the 
evaluation of F(t,t',t) (Ea. (|26[) ). we estimate that the 
integrated noise intensity obtained from Eq. (|27|l entails a 
factor of two theoretical uncertainty for the noise mech- 
anism considered. 

The complex appearance of Eq. 1(27(1 is due to the fact 
that it includes separate terms for the two transitions 
D + *-* D° and D° «-* D~. Such a separation is necessary 
only because the expressions <|13l 114(1 for the activation 
times T + and t_ are slightly different from each other. 
This difference, however, vanishes if one assumes that 
the width of the distributions p(E) and p'(E') is much 
smaller than the width of P(E B ), i.e. 2k B T < AE B . 

The complicated form of the numerator in the right- 
hand side of Eq. I|28|l is the analytic result of the two inte- 
grations over the volume of the "constant resistivity" 
layer. This expression reflects the competition between 
two length scales: the half-thickness of that layer (zt) 
and the screening radius (r s ). In the "thin layer limit" 
zt -C r s , the expression in the square brackets in Eq. lt28|) 

2z 2 

can be approximated by —r- , whereas in the "bulk limit" 
zt r s that expression approaches Azt- 

Now we simplify Ea. ((27|l by taking limits zt <C r s , 
2k B T <C AE B and substituting v(T) = (just an 

approximate numerical fact), which gives 

S v (f) 8^e 4 n D f T (E B )r a s (t (E b )) P{E B )dE B 



V 2 e 2 k B TAEA , 
where 

1 1 



l + 4n 2 pT 2 (E B 



(29) 



1 



ro {Ei 



T+(fJ,,E B )) T-(jI,E B )) 

1.5a, Q (£ B )exp(-:^=#y (:>>()! 



Our actual system is characterized only by the weaker 
inequalities zt < r s , 2k B T < AE B . Therefore, for com- 
parison with experiments we shall still use the original 
formula H27|) . At the same time, for the qualitative anal- 
ysis, we will focus on Ea. (|29[) . but all the conclusions will 
be fully applicable to Ea. (|77|) . 

A remarkable feature of Eq. H29fl . which can be traced 
back to Eq.J7|), is that, even though the resistance noise 
is caused by the fluctuations in the number of conduction 
electrons, the resulting noise spectrum is independent of 
their equilibrium concentration. Furthermore, that spec- 
trum is only weakly dependent on all the defect parame- 
ters involved (see Refii^). In particular, the dependence 



of the prefactor on the concentration of defects tit, is bal- 

— 1/2 

anced by the n D dependence of r S 2 , which then enters 
r 3 s via Eq.lO- 

Like in the BDAE model, the l//-like spectral shape 
generated by formula (|29|) is the result of the broad distri- 
bution of the activation energies P{E B ) 1 but, at the same 
time, formula l|29|l has also two new features, namely: 
(i) the 1/T dependence of the prefactor and (ii) the 
energy-dependent weight r 3 (t(E b )) multiplying P(E B ). 
The second feature is particularly important for extract- 
ing the correct distribution of P(E B ) from experimental 
data. 

For comparison with experiment in Section IVIII 
we will need the integral of S ^P over all fre- 
quencies, which is equal to CV(0). The expres- 
sion for Cy(0) can be obtained by substituting 

(r) 4z T - 3r. s (r) + e~^M(2z T + 3r s (r)) instead of 



1„4 



4z T -3r s (r) + e ~ (2z T + 3r s (r)) 

D( f, t) into Ea. (|27fl . The mathematical structure of this 
expression is similar to Ea. (|27|l but otherwise not very 
illuminating to be written explicitly one more time. In- 
stead, we give an estimate for the noise integral corre- 
sponding to the thin film limit l|29|) . We do it with one 
further simplification. Namely, we replace r s (r(E B )) in 
Ea. (|2~9"|) by the typical value r* = 0.2 /jm and then obtain: 



Syjf) 
V 2 



df = 



2n 2 e A n D rf 
e 2 k B T AE A' 



(31) 



Although the bulk limit zt 3> r s appears to be in- 
compatible with the assumption of independent defect 
fluctuations (to be explained in Section rVIII|l . it is still 
instructive to give an expression for the noise integral in 
this case too. Combining the limit zt ^ r s with all other 
approximations used to derive Ea. ((31|l . we obtain 



f 

■Jo 



Svjf) 
V 2 



df = 



4ir 2 e 4 n D r* 4 
e 2 k B T AE Az T ' 



(32) 



Up to numerical prefactors, both approximations l|31|) 
and (|32|) can be summarized as follows: the noise inten- 
sity is proportional to (j^r^j ■, multiplied by the mean 

squared amplitude of potential fluctuations 6 " e °^ J ^ >s , 
further multiplied by the screening volume of a typical 
fluctuation r* 3 , and, finally, divided by the volume of 
the space, where the defects contributing to the poten- 
tial fluctuations are located. In the thin layer limit that 
volume is roughly 2r s A, whereas in the bulk limit it is 
2ztA. The combination nok B T / AE appearing in the 
estimate of the potential fluctuation should be identi- 
fied with the concentration of "thermally active" defects, 
i.e. those defects that fall in the thermal energy window 
around the chemical potential. 

Now we give various estimates of the prefactor in 
front of the approximate 1// dependence of the spec- 
trum. For the thin film limit, we start from Ea. (|29|l and 



6 



then, assuming r 5 .(r(£ , s)) 
uq(Eb) = uq(Ebo), obtain 



P(E E 



2AE F 



Sv(f) 
V 2 



n 2 e 4 n D rf 



e 2 A AE AE B / ' 
For the bulk limit, the analogous expression is 

Sv(f) 



2Tr 2 e 4 n D r* 



1 



V 2 



e 2 Az T AE AE B f ' 



and 



(33) 



(34) 



Eauation i|34fl admits yet another remarkable simplifi- 
cation, if one substitutes req given by Eq. <|17(l instead of 
r* . (The equilibrium self-screening mechanism described 
by Eg. 1(17) 1 should indeed be a proper description for the 
slower fraction of the defects, which means the lower fre- 
quency part of the spectrum.) In this case 



1 



Sy(f) _ AE 1 

V 2 2Az T n D AE B f ~ N D f 



(35) 



where Njj — IAzttib is the total number of defects in 
the constant resistivity layer. [In the above approxima- 
tion, we assumed AE ~ AEb] Although this result 
looks like a standard statistical factor, its origin is, in 
fact, far more complex. It can be traced back to an acci- 
dental interplay between the resistivity fluctuations and 
the screening mechanism. Given the limitations of the 
"bulk limit" , it is unlikely that formula l|35|) represents a 
clean limit in any realistic system. At the same time, for 
r s ~ zt, this formula still gives a reasonable estimate. 

It is, finally, interesting to observe that, for the num- 
bers characterizing our system, the Hooge's formula™ 
with coefficient a ~ 10~ 3 , would produce 



Sv[f) 



Naf 



an estimate close to that of Eq. (|35|l , (Here Nc is the to- 
tal number of conduction electrons.) This simply reflects 
the fact, that in our case the concentration of conduc- 
tion electrons is about thousand times smaller than the 
concentration of defects. The exact ratio of these con- 
centrations is strongly temperature dependent, as is the 
value of a required to fit the noise intensity in our exper- 
iments. 



VII. COMPARISON WITH EXPERIMENT 

Theoretical spectra computed from Ed. 1(270 are com- 
pared with the experimental ones in Fig. [2Jt (see also 
RefiSi) . The experimental spectra were obtained by sub- 
tracting the zero-current noise from the total noise ob- 
served with V = 50 meV. They were consistent with 
the spectra reported in RefiiS, though the film was 
newly grown. All the theoretical spectra were obtained 
with P(Eb) shown in Fig. \Bp and characterized by 
E B q = V + 0.90 eV and AE B = 0.09 eV. The uncertain- 
ties of the fit for E B o and AE B are 0.05 eV and 0.02 eV, 
respectively. Figure |2t shows r s (Es) obtained by aver- 
aging r s (t + (E, Eb)) and r s (r_(£", E B )) over E and E' 
at T = 434 K. Finally, Fig.JSJi illustrates the role of the 
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FIG. 2: (a) Noise spectra: erratic lines — experiment, smooth 
lines — predictions of Eq. l27ll . The spectra at 340K, 362K 
and 402K are multiplied, respectively, by 30 3 , 30 2 and 30 to 
make them distinguishable. On the real scale all four spec- 
tra fall almost on the top of each other, (b-d) Probability 
distribution P(Eb), the screening radius r a and the product 



P{E B . 
n = 0. 



rf as functions of the barrier energy Eb counted from 



weight by showing P{Eb)t^{Eb)- If P{Eb) were to 
be obtained by fitting the experimental spectra to the 
BDAE modeliS, the result would look like Fig. i.e. 
the maximum would be shifted by about 0.1 eV in com- 
parison with Fig. 

With the above value of AEb , the estimates ()33I34I35|) 
for the prefactor in front of 1// would give respectively: 
9 • 10~ 12 , 3 • lO" 11 and 1 • lO" 11 — all in reasonable 
agreement with experiments. [For the estimates (I34|l and 
((3"5)l we used zt = 0.13 /mi.] 

In Fig. O we present another test of our theory, which 
is independent of the choice of Ebq and AEb- Namely, 
we compare the theoretical and experimental values of 
the integrated noise intensities for the four temperatures 
indicated in Fig. 

The experimental evaluation of the integrated intensity 
of an 1 / /-like noise is a task notorious for its ambiguity. 
Fortunately, in our case, the power law extrapolations of 
all four noise spectra were convergent, which allowed us 
to use the following procedure: 

First, we obtained the lower ends of the error bars by 
integrating the experimental noise spectra only in the 
frequency range of the actual experimental observations. 
Then, the upper ends were obtained by making power 
law extrapolations of the spectra beyond the frequency 
range of observation (up to 10~ 6 Hz for small frequencies 
and 10 8 Hz for large frequencies), and then adding the 
integrals over the extrapolated tails to the lower end val- 
ues of the error bars. Finally, the "experimental" points 
indicated in Fig. |3| were chosen as the middle points of 
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FIG. 3: Integrated noise intensity for the four spectra pre- 
sented in Fig. |21 Empty dots represent experimental val- 
ues, and the solid black dots theoretical values obtained from 
Eg. 1271 . The error bars on the experimental points are ob- 
tained as described in the text 



the above error bars. 

The theoretical points presented in Fig. [3] were ob- 
tained for the full spectrum H27[) as described in Sec- 
tion EH One can also check that the estimate (|31|l gives 
numbers, which are not much different, e.g.: 5.6 • 10~ n 
for T = 340 K, and 4.4 • 1CT 11 for T = 434 K. 

Given the large uncertainties of the extrapolation and 
the theoretical uncertainties indicated earlier, the factor 
of two agreement between the theoretical and experimen- 
tal points in Fig. [3] should be considered as a successful 
consistency test of the noise mechanism proposed. 



VIII. POSSIBLE GENERALIZATIONS 

Now we discuss briefly various modifications required 
in order to generalize our treatment to 1// noise experi- 
ments in other settings involving a-Si:H. 

The application of the present theory to the films with 
coplanar currents^iiiiiSiiii^ and also to the bulk samples 
would encounter an essential complication related to the 
absence of contact layers. In this case, the charge fluctu- 
ations of different defects become correlated via emission 
and capture of the same electron. In such a process, the 
potential fluctuation is simply translated in space from 
one defect to another, which means that, unless the two 
defects have different screening radii, the total resistance 
does not fluctuate at all. Our analysis in Section re- 
vealed one possibility for different defects to have dif- 
ferent screening radii (as a function of their activation 
times). However, even with this possibility, the overall 
effect of the absence of contact layers should be a notice- 
able reduction of the noise intensity. 

The contact layers also play a role in the screening of 
defects. Without them (and without the n-layers) the 
first screening mechanism considered in Section is not 
operational. Therefore, the treatment of screening of the 
faster fluctuating defects becomes more difficult. 



The theory has to be further modified to include the 
hole conduction, if, as a result of doping, the chemical 
potential moves in the middle or below the middle of the 
band gap. An interesting situation may arise, when hole 
resistivity is equal to electron resistivity. In this case one 
should expect a drastic reduction of the noise intensity, 
because the change of electron resistivity induced by po- 
tential fluctuations will be compensated by the change in 
the hole resistivity. 

Analysing experiments, one should also keep in mind 
that creation and annealing of charged defects represents 
an alternative way to potential fluctuations. In partic- 
ular, in a-Si:H such a process may be associated with 
the diffusion of hydrogen atomsi&, which can inhibit or 
expose the dangling bonds. The parameters describing 
this process are not very well known. In our case, the 
mechanism involving the emission and the capture of con- 
duction electrons by already existing defects appears to 
be sufficient, but, in general, creation and annealing of 
defects can lead to a comparable contribution to the in- 
tegrated noise intensity. 

Summarizing this Section, we would like to emphasise 
that all the above complications can still be treated on 
the basis of formula J7J , but different ingredients for the 
evaluation of the potential fluctuations will be required. 



IX. CONCLUSIONS 

In conclusion, we have developed a microscopic theory 
of I// noise in n — i — n sandwich structures of a-Si:H 
and found a very good agreement between this theory and 
our experiments. The noise mechanism proposed should 
be quite general (in particular, in semiconductors): it 
merely requires the presence of charge traps in regions of 
poor screening. The full calculation of the actual noise 
spectrum necessarily involves a fair amount of material 
specific details, which, to some extent, mask the gener- 
ality of our treatment. However, one should remember, 
that our calculation is based on formula Q, which ap- 
pears to have quite a broad range of applicability. 

From a different perspective, our analysis introduces a 
new method of characterizing the deep defects in a-Si:H. 
In particular, our finding, that Ebo — E c q = 0.27 eV, 
which is substantially greater than 0.1 eV - - the 
scale of potential fluctuations seen in the drift mobil- 
ity experiments^, may indicate that the defects are not 
distributed randomly in the background of the medium- 
range random potential but, instead, located in the re- 
gions of greater local strain, i.e. around the peaks of that 
potential. 

Our result also limits the spectral intensities due to 
all other possible noise mechanisms to the difference be- 
tween the experimental spectra and the spectra given 
by Ea. (|27[l . Even with the freedom of varying Ebo and 
AEb, this is a very strong constraint. 
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